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K-SHELL  YIELD  SCALING  LAW  FOR  CONVENTIONAL  PRS  LOADS 


I.  INTRODUCTION 

An  interesting  figure  .of  merit  often  .used  to  characterize  the 
performance  of  Z-pinch  devices  such  as  BLACKJACK  5  at  Maxwell  Laboratories, 
Inc.,  and  PITHON,  of  Physics  International  Corp.  is  the  energy  released  in 
K-shell  radiation  during  a  discharge  with  a  given  element  as  the  plasma 
load.  On  both  PITHON  and  BLACKJACK  5,  K-shell  yields  have  been 
experimentally  measured  for  load  elements  ranging  from  neon  (Z=*10)  to 

krypton  (Z=36).  To  date,  BLACKJACK  5,  and  DOUBLE  EAGLE  with  a  maximum 
power  of  approximately  8-10  TW,  represent  the  leading  edge  in  the 
technology  of  such  devices. 

For  theoretical  analyses  of  the  effects  produced  by  these  generators, 
a  variety  of  tools  and  strategies  have  been  developed,  ranging  from  simple 
"zero-dimensional"  implosion  models1  to  detailed  MRHD  (Magneto  Radiation 
HydroDynaraic)  codes  containing  atomic  ionization  dynamics  and  radiation 
transport  capabilities.^’^  While  the  latter  detailed  models  are 

indispensible  toward  achieving  an  understanding  of  the  complicated  physics 
issues,  it  is  oftentimes  appropriate  to  appeal  to  simplified  models  in 

order  to  gain  physical  insight  and  order  of  magnitude  estimates  of  the 
various  effects  in  a  noninteractive  mode.  In  this  Memo  Report  it  is  our 
intent  to  bridge  the  gap  between  the  overly  complicated  and  the  overly 
simplistic  approaches  by  considering  a  simple  dynamic  model  for  the 

implosion  coupled  with  atomic  physics  and  empirical  data  based  on 

experimental  observation.  Our  principal  result  is  the  establishment  of  a 
predictive  capability  which  will  serve  as  a  guide  for  estimating  generator 
power  requirements  for  substantial  K-shell  radiation  for  conventional  PRS 
loads . 

II.  BASIC  MODEL  ELEMENTS 

1  —  ? 

The  equation  of  motion1  J  describing  the  implosion  of  a  plasma  annulus 
when  the  magnetic  pressure  greatly  exceeds  the  plasma  pressure  is 

2  -7  2 

d  r  -10  I(tT 

m  — r  -  - - -  ( 1 ) 

dt 
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where  m  is  the  annulus  mass  per  unit  length  (kg  m”*),  r  is  the  radius  (m) , 
and  I  ( t )  is  the  time-varying  current  (amps).  For  a  wide  variety  of 
machines  from  GAMBLE-II  to  BLACKJACK  5,  the  main  current  pulse  is 
adequately  described  by  a  linear  ramp  I(t)  ■  -IQt  where  IQ  Is  the  current 
rise  slope  in  amps  sec”*.  Fig.  1  shows  the  actual  current  pulse  for 
BLACKJACK  5  shot  177  and  the  corresponding  excellent  linear  fit  to  I(t) 
during  the  main  run-in  phase.  Equation  (1)  with  the  linear  current  fit 
becomes 


m 


-7  2  2 

-10  I  t 
o 


(2) 


Of  course,  the  resisting  plasma  pressure  late  in  the  run-in  will 
render  this  simple  picture  inaccurate.  However,  considerable  evidence  from 
various  machines  strongly  indicates  that  Eq.  (2)  obtains  until  the 

plasma  has  collapsed  to  about  10%  of  its  initial  radius,  at  which  time  the 
implosion  abruptly  ceases  and  the  main  radiation  pulse  commences.  Our 
description  is  based  upon  this  empirical  observation.  Eq.  (2)  may  be 
universalized  by  adopting  the  following  dimensionless  variables,  where  RQ  * 
initial  radius  in  meters. 


x  =  r/R. 


q  =  1.78  x  10“2  (I  /R  )1/2  m“1/4  t 
o  o 


(3a) 

(3b) 


Eqs.  (3a)  and  (3b)  combined  with  Eq.  (2)  lead  to 
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d“x 


o 

q" 


(4) 


The  boundary  conditions  applied  to  Eq.  (4)  are:  at  t=0,  x=l,  and  all  q- 
derivatives  of  x  are  0,  since  l(t)=0  at  t=0.  To  our  knowledge,  no  analytic 
solution  of  Eq.  (4)  exists,  but  a  power  series  solution  is  obtainable. 
This  series  solution  converges  too  slowly  to  be  of  real  practical  value 


x(q) 


4  8  12 

_  3 _ 3 _ _3 _ 

12  672  15,650 


(5) 
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Fig.  1.  Current  trace  for  BLACKJACK  5  shot  177,  along  with  the  linear  fit 
adopted  to  describe  the  major  run-in  phase.  Current  trace  provided 
courtesy  of  Maxwell  Laboratories,  Inc. 
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The  exact  numerical  solution  of  Eq.  (4)  is  presented  in  Fig.  2.  At 
our  "magic"  r/R0  ■  x  ■  0. 1,  where  the  implosion  is  deemed  to  cease  and 
thermalization  begins,  the  time  variable  q  -  1.70,  and  the  implosion 
velocity  dx/dq  ■  3.3.  Key  physical  variables  at  this  point  are  the  kinetic 
energy  per  ion  and  the  implosion  time.  Since  q  «  1.70  the  implosion  time 
(sec)  is  easily  obtained  from  Eq.  (3b)  as 

R  1/2 

t.  =  95.6  (-=2-)  m^4  sec.  (6) 

imp  V 

To  obtain  the  kinetic  energy  per  ion  at  x  ■  0.1,  note  the  atomic  weight  A  as 
1.32  Z1*15  where  Z  is  the  element's  atomic  number.  Straightforward 
algebraic  manipulation  of  the  variables  leads  to 

2 

1. 64xl0~12  R  I  A  f^-1 

(ke)  .  - - 2_° - eV  (7a) 

lon  /m 

or,  since  dx/dq  =  3.3  at  x  *  0.1,  and  A<S,1.32Z**^ 

2. 36xl0_l 1  R  I  Z1,15 

(KE)  -  - 2-2 -  eV  .  (7b) 

Ion  , 

/m 

2  3 

Conventional  loads  are  usually  designed  >  so  that  the  implosion  time, 
Eq.  (6),  matches  the  time  of  peak  current,  I_.  In  some  cases  the  match  is 
only  approximate,  but  even  a  10-20%  deviation  does  not  substantially  change 
our  results.  The  current  slope  IQ  =  lp/timp’  -sing  Eq.  (6)  for  the 
implosion  time,  t^mp,  leads  to  a  relationship  between  I0  and  Ip 

1.09xl0~4  I  2 

I  -  - 2—  amps  sec  .  (8) 

°  R  /m 
o 

Substituting  for  IQ  in  Eq.  (7b)  we  finally  obtain  the  following  expression 
for  kinetic  energy  per  ion  at  10%  of  the  initial  radius,  where  we  assume 
that  thermalization  and  conversion  to  radiation  of  the  kinetic  energy 
occurs,  viz. 
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eV 


(9) 


2.5xl0~15  I  2  Z1,15 


(KE) 


ion 
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III.  BASIC  ATOMIC  PHYSICS  ENERGETICS 

If  the  K-shell  of  a  given  element  is  to  be  excited,  an  expended 
minimum  of  energy  which  must  be  imparted  to  each  imploding  ion  is  that 
needed  to  strip  all  but  one  or  two  electrons  from  the  ion  plus  the  energy 
needed  to  maintain  each  stripped  electron  at  an  appropriate  temperature 
conducive  to  K-shell  dominance.  Of  course  inefficiencies  and  losses- 
particularly  radiative  -  imply  that  in  general  much  more  energy  than  this 
is  needed.  However,  much  insight  into  basic  Z-scaling  may  be  obtained  by 
consideration  of  these  basic  quantities. 

The  temperature  at  which  copious  K-shell  radiation  has  been  observed 
in  Z-pinches  is  given  approximately  by 

z  2 

T  -  1300  (y^j  eV  .  (10) 

We  approximate  the  number  of  stripped  electrons  as  equalling  Z.  Therefore, 
to  maintain  the  stripped  electrons  from  an  ion  of  atomic  number  Z  at  this 
temperature  requires, 


E  ■  4  Z2  eV/ion.  (11) 

m 

The  energy  needed  to  strip  the  electrons  from  each  ion  is  obviously  also 
important  and  non-negligible.  We  shall  assume  that  the  plasma  is  50% 
hydrogen-like  and  50%  helium-like.  The  ionization  potentials  for  each 
element  are  available  from  many  references;  we  have  used  the  data  of 
Carlson  et  al.®  The  energy  needed  to  strip  an  ion  of  atomic  number  Z  to  a 
mean  charge  state  of  Z— 1.5  is  given  by  the  following  fit  -  which  for  the 
approximate  treatment  of  this  report  has  negligible  error. 

2  6  A 

E  -  3.83  Z  *  eV/ion  (12) 

s 
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Obviously,  the  expended  minimum  of  energy  per  ion  required  is  given  by 

Em+Eg.  We  now  inject  another  empirical  observation.  A  critical  question 

is,  in  practice,  how  much  energy  per  ion  is  required  for  copious  K-shell 

radiation  in  Z-pinches?  BLACKJACK  5  shot  177,  an  aluminum  wire  array 

implosion  which  produced  -  20  kJ  of  K-line  radiation  was  one  of  the  best 

diagnosed  shots  in  the  BLACKJACK  5  series.  From  the  current  trace, 

spectroscopic  diagnosis,  and  known  characteristics  of  the  load,  it  was 

determined  that  3.6  x  10^  eV  was  needed,  on  average,  to  produce  a  K-shell 

ion  in  this  shot.  This  number  is  30  times  E_  +  E_.  We  now  make  the 

m  s 

assumption  that  for  all  Z,  30  times  (Em+Eg)  will  be  required  for  efficient 
K-shell  production.  Such  a  large  factor  is  plausible  when  one  considers 
the  numerous  sinks  for  generator  energy  other  than  production  of  K-shell 
radiation.  Such  processes  include  joule  heating  early  in  the  implosion, 
radiative  losses  other  than  K-shell,  thermal  conduction  losses,  and 

inductive  losses.  This  is  not  unreasonable  as  a  first  approximation  since 
it  allows  radiative  and  other  losses  to  scale  upward  with  Z  in  the  same 
fashion  as  the  kinetic  and  ionization  energies,  which  is  intuitively 

reasonable.  This  factor  of  30  can  be  viewed  optimistically,  since  it 

indicates  much  room  for  improvements  in  efficiency  of  these  devices.  This 
assumption  means  that 

(KE).  =  30  (4Z3  +  3.83  Z2'64)  as  187  Z2’95  eV/ion  (13) 

ion 

where,  in  Eq.  (13)  the  additional  simplification  of  the  power-law  sum  to  a 
single  power  law  does  not  result  in  substantial  error  for  present 

purposes.  Perhaps  the  most  important  single  step  in  our  development  now 
follows,  namely  equating  Eq.  (13),  the  required  kinetic  energy  per  ion,  to 
Eq.  (9)  the  actually  imparted  kinetic  energy  per  ion 


2.5xlO~15  I  2Z1,15 


m 


=  187  Z 


2.95 


which  yields 


m  =*  1.34x10  I  2  Z  kg/m 
P 


(14) 


(15) 
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for  the  mass  per  unit  length  which  can  be  stripped  to  the  K-shell  for  an 
element  of  atomic  number  Z  with  a  peak  machine  current  I  amps. 

An  ion  density  may  be  obtained  from  Eq.  (15)  if  the  radius  of 
implosion  is  known.  No  systematic  variation  of  the  pinch  diameter  with  Z 
has  been  seen;  we  adopt  1.34  mm  as  an  approximate  average  which  was 
observed  on  successful  shots  using  pinhole  imagery,  and  assume  future 
machines  will  produce  similar  pinches.  Given  the  total  mass,  the  mass  per 
ion  (2.19xl0~^7  Zl-15  kg),  and  the  cylinder  radius  the  K-shell  ion  density 
is  obtained  easily  as 

N  -  4.33xl09  I  2  z“2*95  cm-3  .  (16) 

Ip 

From  extensive  spectroscopic  density  and  temperature  diagnostics,  we 

have  determined  that  the  yield  on  previous  shots  obeys  the  following 

—3 

function  of  temperature  T(eV)  and  ion  density  Nj(cm  ). 


Y(kJ)  a  142(£rJ  [' 


,  _  , -19/ Z  .2 
3.7x10  (yg-1 


■l1*  9  exp  [  -2xl03  fyjr-J2  /t] 


.  min  [1.0,  2xl03  \j^y  /t]  . 


We  emphasize  that  Eq.  (17)  is  not  based  on  any  hydrodynamics;  only  on 
spectroscopic  measurements  of  Nj  and  T  coupled  with  knowledge  of  the  yield 
on  the  diagnosed  shots.  The  1.9  power  law  of  Nj  (instead  of  the  usual  2.0) 
represents  an  approximate  correction  for  opacity  effects.  Letting  T  =  1300 
(Z/18)^  eV  (Eq.  10),  Eq.  (17)  is  greatly  simplified  to 

Y  =*  3.66  x  10“34  N**9  Z”1*8  KJ  .  (18) 

Substituting  for  Nj  from  the  zero  dimensional  hydrodynamics  -  atomic 
energetics  model,  given  by  Eq.  (16)  we  obtain 

Y(kJ)  =■  7.46  x  10-16  I3*8  Z-7*405  KJ.  (19) 

P 


For  BLACKJACK  5  the  power  in  TW  (P)  is  given  approximately  as  3.86  x 
—  13  2 

10  I _  ,  resulting  in  the  following  alternative  form  of  Eq.  (19). 


KJ 


(20) 


Y(KJ)  =  2.87  x  108  P1,9  Z-7*405 


IV.  CONCLUSIONS 


We  have  combined  scaling  laws  for  empirically  observed  yield  as  a 
function  of  spectroscopically  diagnosed  temperature  and  density  with  basic 
atomic  physics  energetics  and  a  standard  zero-dimensional  equation  of 
motion  for  an  imploding  Z-pinch  load  to  obtain  a  prediction  for  yield  for 
any  element  as  a  function  of  peak  current  (Eq.  19)  or  of  BLACKJACK  5  -  type 
machine  power  (Eq.  20).  The  predictions  of  Eq.  20  are  listed  in  the 
accompanying  table,  along  with  the  actual  BLACKJACK  5  optimized  yields. 
Clearly  this  model  does  a  very  good  job  of  describing  previous  results. 
When  it  is  applied  to  more  powerful  machines  contemplated  for  future 
construction  several  important  implications  are  seen.  First,  a  25-TW 
generator  would  render  a  conventional  argon  load  a  very  prolific  PRS,  with 
K-shell  yields  probably  exceeding  50  KJ.  Somewhat  larger  powers  in  the  40- 
50  TW  range  would  be  needed  for  similar  success  with  titanium.  Finally, 
conventional  loads  consisting  of  iron  or  krypton  will  require,  for 
substantial  K-shell  line  radiation,  superclass  generators  in  the  50-100  TW 
range. 

Further  review  of  the  table  shows  that  by  doubling  the  generator  power 
the  K-shell  radiative  yield  quadruples.  This  suggests  that  more  current 
would  couple  into  the  load  to  heat  the  plasma  and  that  the  plasma  becomes 
hotter,  i.e.,  the  excitation  temperature  associated  with  the  K-shell 
emitting  region  is  higher.  The  question  of  whether  more  current  actually 
will  couple  to  the  load  is  somewhat  speculative  in  the  sense  that  there  are 
no  direct  measurements  of  current  in  the  vicinity  of  the  load  but  only 
measurements  of  the  return  current.  .Also  evident  from  the  table  and  Eq. 
(20)  is  that  there  is  a  tremendous  price  to  pay  in  terms  of  energy 
expenditure  in  stripping  moderate  and  high  Z-elements  down  to  and 
maintaining  the  K-shell  plasma.  The  sequence  of  events  leading  to  K-shell 
emission  involves  both  heating  and  compressing  the  plasma.  Ifeture  is  not 
kind  to  these  schemes.  It  might  be  more  efficient  to  devise  schemes  and 
design  loads  that  first  compress  and  then  heat  the  plasma. 


Ic  Is  certainly  possible  that  advanced  concepts  for  improving  both 
machine  technology  and  the  pinch  itself  through  better  understanding  of  the 
influence  of  current  pulse  duration  and  crucial  plasma  instabilities  may 
favorably  alter  these  predictions.  For  instance,  a  number  of  novel 
experiments  involving  pulse  shaping  through  such  devices  as  the  plasma 
erosion  opening  switch  are  being  investigated  to  determine  whether  it  is 
possible  to  enhance  K-shell  yields  in  conventional  PRS  loads.  Unless  a 
greater  bulk  of  the  plasma  can  be  heated  and/or  further  compressed  through 
pulse  shaping  it  is  difficult  to  conceive  how  the  energy  expenditure  for 
stripping  and  maintaining  a  K-shell  plasma  can  be  reduced  from  the 
efficiency  factors  discussed  here. 

TABLE:  Prediction  of  Eq.  (20)  for  K-shell  yield  in  kilojoules  as  a 

function  of  atomic  number  Z  and  generator  power  in  TO.  A  BLACKJACK  5-type 
machine  with  a  load  length  of  3  cm  is  assumed.  The  asterisk  indicates 
where  the  prediction  is  probably  invalid  due  to  burnthrough  of  the  K- 
shell.  The  numbers  in  parentheses  are  Maxwell  Laboratories'  optimized 
results,  courtesy  of  M.  Gersten. 
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* 
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